De-Rotation of the Cosmic Microwave Background Polarization: Full-Sky Formalism 
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Mechanisms have been proposed that might rotate the linear polarization of the cosmic microwave 
background (CMB) as it propagates from the surface of last scatter. In the simplest scenario, the 
rotation will be uniform across the sky, but the rotation angle may also vary across the sky. We 
develop in detail the complete set of full-sky quadratic estimators for the rotation of the CMB 
polarization that can be constructed from the CMB temperature and polarization. We derive the 
variance with which these estimators can be measured and show that these variances reduce to the 
simpler flat-sky expressions in the appropriate limit. We evaluate the variances numerically. While 
the flat-sky formalism may be suitable if the rotation angle arises as a realization of a random field, 
the full-sky formalism will be required to search for rotations that vary slowly across the sky as well 
as for models in which the angular power spectrum for the rotation angle peaks at large angles. 

PACS numbers: 98.80.-k 
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I. INTRODUCTION 

Great strides have been made during the past decade 
in obtaining precise maps of the cosmic microwave back- 
ground (CMB) temperature and polarization \w. Yet 
there are still increasingly precise data to come Q. In- 
flation 0, which has passed a number of tests so far, will 
be tested further 0] , and we may begin to understand the 
physics responsible for inflation. Moreover, there will be 
a number of different types of departures from the stan- 
dard scenario that can be detected with new measure- 
ments. 

The CMB polarization pattern can be decomposed into 
a gradient component (E modes) and a curl component 
(B modes) HQ. A given inflationary model makes spe- 
cific predictions for the power spectra of these compo- 
nents. For example, if the model predicts no primor- 
dial gravitational waves, then there is no B mode at the 
CMB surface of last scatter. If there is a primordial 
gravitational-wave background with a sufficiently large 
amplitude, there will be B modes at the surface of last 
scatter at a detectable level Q , and with a characteristic 
power spectrum. 

Since the post-recombination Universe is transparent 
to the CMB, it is natural to assume that the polariza- 
tion pattern that we see is the polarization pattern at the 
surface of last scatter. However, exotic mechanisms may 
rotate the linear polarization as it propagates through the 
Universe. For example, if cosmic acceleration is due to a 
quintessence field 4>, and if that field couples to the pseu- 
doscalar FF of electromagnetism, then the time variation 
of 4> leads to a rotation of the linear polarization by an 
angle a [f|. If is spatially homogeneous, then the ro- 
tation angle a is uniform across the sky. In this case, 
there is a characteristic parity-violating EB correlation 
induced in the CMB polarization Q . Null measurements 
of this effect constrain a to be smaller than a few degrees 

It may well be, however, that <j) has spatial variations, 



in which case the rotation angle a (ft) will vary with po- 
sition ft on the sky [l2j . There may also be mechanisms 
involving coupling of photons to dark matter that cause 
a spatially- varying a(fi) [l3[. In either case, the EB 
power spectrum would include the dependence on a (ft), 
which would differentiate this scenario from the uniform- 
rotation case. For example, the EB power spectrum even 
vanishes fl2j if a has a dipolar variation across the sky, 
in which case the positive EB correlation on one half of 
the sky is cancelled by a negative EB correlation of equal 
magnitude on the other half. 

Ref. [3| showed that the rotation angle a (ft) can be 
reconstructed, as a function of position on the sky, from 
the measured polarization map. If primordial perturba- 
tions are assumed to be Gaussian (as predicted by infla- 
tion), then the rotation introduces a characteristic non- 
Gaussian signal in the polarization map. Measurement 
of this non-Gaussian signal then provides the rotation 
angle. More specifically, if primordial perturbations are 
Gaussian, then the spherical-harmonic coefficients Ei m 
and Bi m of the polarization map are statistically inde- 
pendent for different I and m. However, the rotation 
introduces off-diagonal correlations, i.e. correlations be- 
tween different Im pairs. 

Ref. [IH re-visited this proposal. The authors worked 
out the off-diagonal correlations that would be induced 
in the flat-sky limit, and evaluated numerically how well 
the small-angle analogues of the spherical-harmonic coef- 
ficients (Xlm of the rotation angle could be measured with 
future satellite experiments like Planck and CMBPol 

In this paper, we work out in more detail the full-sky 
formalism of Ref. [14|. We write down the minimum- 
variance estimator c\lm for the rotation-angle spherical- 
harmonic coefficients that can be obtained from a full-sky 
CMB temperature/polarization map. We derive expres- 
sions for the variance with which the aru can be deter- 
mined and show that they recover, in the large-L limit, 
the correct flat-sky expressions. While a flat-sky anal- 
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ysis may be suitable if the rotation-angle power spec- 
trum peaks at small scales, the full-sky formalism will 
be required to maximize the sensitivity in models, such 
as that in Ref. [T3|, where the signal- to- noise peaks at 
low L (see, e.g., Fig. 1 in Ref. [15]). There has also 
been growing attention recently to the possibility that 
there may be variations in fundamental fields over dis- 
tance scales comparable to, or larger than, the horizon 
(perhaps remnants of the pre-inflationary Universe) [l8| . 
Observationally, these entail searches for departures from 
homogeneity/isotropy or departures from statistical ho- 
mogeneity/isotropy in the CMB or large-scale structure 
[l9l . |20| . The full-sky formalism we present here can be 
used to search for the low-L (e.g., L = 1, 2, 3, • • • ) mo- 
ments of a(n) that may arise if <f> has long-wavelength 
fluctuations, in addition to, or instead of, the higher-!/ 
modes that can also be probed with a survey of a smaller 
region of the sky 1 . 

The plan of the paper is as follows. In §11, we derive 
the temperature/polarization correlations induced by ro- 
tation. In §111, we construct the complete set of quadratic 
estimators for the rotation of the CMB polarization that 
can be constructed from the CMB temperature and po- 
larization maps. We also determine the variances with 
which these estimators can be measured with a full-sky 
CMB map. §IV derives the flat-sky limit of our variance 
expressions and compares them with previous work. §V 
presents numerical results for the variances. Conclud- 
ing remarks are provided in §VI. Appendix A provides 
some useful formulas. Appendix B shows how the effects 
of weak lensing and rotation can be distinguished geo- 
metrically. Appendix C shows that our variance expres- 
sions retrieve successfully the expressions for a uniform 
rotation angle. And Appendix D illustrates the relative 
contribution of different IV pairs to the rotation-angle 
estimator. 



II. CORRELATIONS INDUCED BY ROTATION 

In this Section, we derive the temperature/polarization 
correlations that are induced by a post-recombination ro- 
tation of the polarization. 



A. Induced Modes 

We begin by recalling that the Stokes parameters Q(n) 
and U(n), as a function of position n on the sky, are 



components of a symmetric trace-free 2x2 tensor, 2 



P ab (n) 



Q(n) -U(n)sini 
2 V — U(n) sinf — Q(n)sin 2 



(1) 



This tensor can be expanded in terms of tensor spherical 
harmonics Y,¥ n . (n) and Y,f , , (n) in the usual fashion 

SI, 



oo X 



P ab (n)=J2 E i E lrnY ( f m)ab (n) + B lm Y ( f m)ab (n)}, 

1=2 m=-l 

(2) 

where Ei m arc the E-modc tcnsor-sphcrical-harmonic 
coefficients, and Bi m are the B-mode tensor-spherical- 
harmonic coefficients. 

We now suppose that the polarization at the surface of 
last scatter is a pure E mode; we comment on the validity 
of this assumption below. Then, a small rotation a(h) 
induces a change to the polarization, 



SP ab (n) = 2a(n)P; b (n), 



where UA 



oo I 



(3) 



(4) 



1=2 m=-l 



A pure E mode thus gets rotated into a B mode. Note 
that Eq. ([4| is valid only in the limit of small rotations, 
a(n) <C 1. Given that, the B mode will be small com- 
pared to the E mode, which satisfies existing empirical 
constraints. 

Rotation is a scalar field on the sky and can thus be 
expanded in terms of spherical harmonics, 



(5) 



LAI 



where ulm arc the coefficients of the expansion. 

Using Eqs. d3j) , (J4j) , and ([5]), the B mode induced by a 
small rotation angle from a pure E mode can be expressed 
as 



8Bi m = J d?iSP ab (h)Y ( f^ ab (h) 
= 2 E E a LMEi 2m2 / dhY^ b Y {LM) Y { 

r i/ i *f 



B 



(6) 



LAI (2Tfl2 

Similarly, the induced E mode is 

rE*,ab/ 



SE im = / dn5P ab (h)Y^ a (n) 



2 EE a «' £, « / dfiY S™) abY ( LM ) Y { 



B 

{l2Tn2)ab ' 



(7) 



LAI l 2 m 



2 Note that the definition of P a i, here differs from that in Refs. || 
1 The full-sky formalism is exact; note, however, that it is not so that the normalization of the power spectra agree with those 

computationally more demanding than the flat-sky calculations. of Refs. 
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We further develop the last two expressions using Q 

Y (fm)ab = -ir(X{lm):ac£ C b + Y( lm ). bc e C a) , (8) 



Y {fm)ab — N l(Y(lm):ab ~ -^9abY(l m ) :c C ) , (9) 

where g and e are, respectively, the metric tensor and the 
Levi-Civita tensor on a unit 2-sphere (see Appendix A); 
a colon denotes a covariant derivative on the 2-sphcrc; 
and N\ is given by 



N, 



'2(1-2) 
"0 + 2)1 



(10) 



From Ref. [U, we can express the double derivatives in 
terms of spin-2 spherical harmonics ±2Yn m ) (fi) as 



1(1+1) ^ 

r (Im):at — ^ * (lm)9ab 



1 (J + 2)! 



2 ' 2 y (Z — 2)! (11) 

x [2Y(i m ){m + ®m + ) +_ 2 F (iTO) (m_ ®m_)] a &, 

where 



1 



m± = -p(ee + ie^). 
V2 



(12) 



Combining Eqs. (JTTJ and (O with Eqs. ©, ©, and 
(|10[) , we can express the tensor spherical harmonics Y, 



in terms of spin-2 spherical harmonics as 3 



E,B 
(lm) 



Y, 



(lm) — 

^( +2 Y-- 2 Y) ^sm(9)UY+ +2 Y) \ 
^sm(9)UY+ +2 Y) ^^{9)UY- +2 Y) ' 



(13) 



Y 



(lm) 



% +2 Y+_ 2 Y) i^ S m(6)(- 2 Y- +2 Y) 
^ sin(0)(_ 2 y- +2 Y) sin 2 (0)(_ 2 y+ +2 Y) 



Using this result, we obtain 

yB*,abyB 
1 (lm) I (l 2 m 2 )ab 



1 , , (15) 

l}(-2Y* lm ) x_ 2 Y(i 2m2 ) + +2 Y(; m) x +2 Y(i 2m2 )), 



and 



yE*,ab-yB 
(lm) I (l 2 m 2 )ab 

(16) 

- -^{+ 2 Y* lm) x +2 Y(i 2Tn2 ) -_ 2 y"( im ) x_ 2 Y(7 2m2 )). 



We suppress the (Zm) indices for Y in this formula. 



The next step is to use Eqs. (fl5j) and (| 16)) in order 
to rewrite the integrals in Eqs. © and ([7]) in terms 
of Wigner 3j symbols. We use Ref. [HI], noting that 
the spin-spherical harmonics of zero spin are the regu- 
lar spherical harmonics, oYa m ) = Ya m \. Let us first look 
at the case of the induced B mode, where we have 



dhY,?*l ab Y (r ,M,Y? 



(lm) I (LM)^(l 2m2 )ab 



I L l 2 

-2 2 



We now define 4 



cLM _ c_iy 

^lml 2 m 2 V / 



-in 

1 L l 2 

2 0-2 



I L h 
-to M m 2 



(2/ + l)(2L + l)(2i a + l) 



47T 



and 



l L l 2 
-m M m 2 



1 L h 

2 0-2 



(17) 



(18) 



(19) 



Also, note that, due to the properties of the Wigner 3j 
symbols, 5 the sum in Eq. p7|) vanishes, unless l + l 2 + L = 
even. Replacing Eqs. l[T7]). ([T8]> . and HU) in Eq. ©, we 
come to a relatively simple expression for the rotation- 
induced B mode, 



5B bm - 2^ ^2 <XLMEi ama €t ml . 

LM lo,m<z 



-LM 



(20) 



where the only non-zero terms in the sum are those that 
satisfy I + l 2 + L = even. 

Similarly, for the case of the rotation-induced E mode, 
again using the properties of the Wigner 3j symbols, 6 we 
have 



dnY^ b Y (LM) Y ( l m2)ah = 



ml 2 m 2 ll 2 



(14) and then replacing this in Eq.([7]), we get 



6E lm - 2«^ 2J a LAlEl 2 m 2 £,lm X l 2 m 2 H < 



(21) 



(22) 



LAI l 2 m 2 



Note that in Eqs. ([2"T]) and (|2"2")l the only non-zero terms 
are those that satisfy I + l 2 + L = odd. 



4 Note that the definitions of , and H}i, differ from those in 
Ref. I14H . We define these quantities in this way to avoid division 
by zero. 

5 See Appendix A: changing the sign on all three m's, brings up a 
factor of (-l) l + l i+ L . 

6 See Appendix A: if the sum of m's does not vanish, the value of 
the symbol is zero. 
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B. Induced Correlations 

Using the induced modes, derived in the previous 
subsection, we derive four correlations that are modi- 
fied/induced by the rotation: EB, EE, TB, and TE (there 
is also a BB correlation, but it is higher order in a, and 
thus neglected). Every mode that we detect will contain 
the sum of the primordial (at the surface of last scatter) 
and the rotation-induced component, 



Ehn — Elm.Q + SEl m , 
Bl rn = SBl rn . 



(23) 



We use the definitions of the EE and TE power spectra, 

(Elm,oE*t m i ) = Cf E 'Sll'Smm' , 
(Elm,oTi/ m r Q ) = Cf E 'Sll'S mm ' . 



(24) 



Finally, we obtain the expressions for the four correla- 
tors, to first order in a. Using Eqs. (HOI) and l|2"2" j) -([2"g ]) . 
we obtain 



1 



LAI 



=a 00 Cf E SwS 



(25) 



+2E 



L 



r*EE c LM ttL 



L>1 M=-L 



and the rest are 

(E lm E;, m ,) = c EE s w s, 



LM ttL 

Iml'm' IV 1 



(26) 



LM 



{BlmTi*, m> ) - 2 U L MCj E i{ml' m 'HlV i ( 27 ) 



LAI 



(EimT^mi) — Cj E 1 5ll'8 mm ' 

Iml'm' -^W ' 



LM 



(28) 



Note that for EB and TB, the sum is taken over the 
terms that satisfy I + I' + L = even and in EE and TE 
over I + I' + L = odd. 



A 


Z ll> 


BE 


°i' 


EB 




EE 


-uc EE - cF E ) 


BT 




TB 


cr 


ET 


-icF 


TE 





E3, for the 



TABLE I: The quantities Zf v , defined in Eq. 
various modes A. 



show however, rotation induces off-diagonal correlations; 
i.e., correlations between spherical- harmonic coefficients 
of different lm and I'm'. 

While the correlations of specific lm-l'm' pairs depend 
on the azimuthal quantum numbers m and m' , they can 
be parametrized in terms of the rotational invariants (i.e., 
independent of m) , 



D iv M ' = ^ a LAiZa,H E v 



L 

IV ') 



(29) 

where the quantities Zjf, are given in Table 1; they are 
obtained from Eqs. (EHJ)- The EB correlator, for a 
given IV pair with I ^ I', is different than the BE cor- 
relator. We thus consider both BE and EB (and sim- 
ilarly for TE/ET and TB/BT) and then restrict our 
sums to I' > I to avoid double-counting. We thus have 
for I > V, A = {BE,EB,TB,BT,TE,ET,EE}, while 
for I = I', we have A = {BE,TE,TB, EE}. With 
these shorthands, the part of the XX' correlators (for 
{X, X'} = {T, B, E}) can be written as 



{Xi m {X Vm ,) ) = 2_^ u iv 6 



Iml'm' 



(30) 



LM 



We now suppose that we have spherical-harmonic co- 
efficients T™ ap , £™ ap , £™ ap , obtained from a CMB 
temperature/polarization map. These receive contribu- 
tions from the true signal on the sky, reduced by the 
Z-space window function Wi = exp(— Z 2 cr 2 /2), where 
°~b = ^fwhm/v81n2 = 0.00741 (0f w hm/l°) for a Gaussian 
beam of width #f w hm, and a contribution from detector 
noise. The predictions for the rotational invariants for 



D LM,XX' Wi _ p oUowing 



the map are U u , 

Rcfs. [l4|, EH, EH, the minimum- variance estimator for 

i T-.LM. XX' .man . 7 

each D u , F is 



D 

u iv 



LM,XX',map 



1 r-iL \ — 1 \ r ymapy/,map,*>IM 



III. ESTIMATORS FOR THE 
ROTATION- ANGLE SPHERICAL-HARMONIC 
COEFFICIENTS 

If we assume that the primordial CMB tempera- 
ture/polarization pattern is a realization of a statisti- 
cally isotropic Gaussian random field, then the spherical- 
harmonic coefficients (T; m , Bi m , or Ei m ) for the pri- 
mordial field are all uncorrelated. As Eqs. (f2"5")) - (f2"5|) 



where we have used (see Appendix A), 



t LM 
■slml'm' 



in 



(31) 



2 = G\ v = (g[±lg±ll. (32) 



7 Note that the definition of Gh, differs from that in Refs. [I3.H9||. 
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Recall also that for EB and TB, only the terms in Eq. (|3T|) 
that satisfy I + I' + L =even are non-vanishing, while for 
EE and TE only l + l' + L =odd terms are non- vanishing. 

The variances with which each ^ ' map can be 

measured can also be calculated. Moreover, measure- 
ments of different Df v ' ' map will be correlated, for 
the same IV. We write the covariances between the 
L>f l f /I,XX ,map in terms of the quantities, 



nil /~iL 

U AA' = U W 



D 



L M , A,map j^L M , A' ,map \ 



£yLM, A, map \ 



IM,A'map\ 



^iv 



(33) 



In principle, C\ A t is a 7 x 7 matrix (in the AA' space) 
for I 7^ V and 4 x 4 for I = V. However, the matrix is 
sparsely populated in the AA' space; it can be written 
in block-diagonal form, since the EB and TB correlators 
arc non- vanishing only for l + l' + L =evcn while those for 
TE and EE are non-vanishing only for I + I' + L =odd. 
We present explicit expressions for the relevant entries of 
C\ A , below. 



A. Minimum- Variance estimators for the 
rotation-angle coefficients 

We now write down the minimum-variance quadratic 
estimator that can be constructed for the rotation- 
angle coefficients a^M from the measured tempera- 
ture/polarization map. 

To begin, we recall that each D^f 4 ' ,map (i.e., each IV 
and AA' , for a given LM) provides a measurement of 

ollm, through (o.lm)u, 
the quantities 



_ j\LM, A, map 
U ll 



/Fj^; A , in terms of 



p L.A 



2Z$H l L l ,W l W l < 



(34) 



The variances and covariances of the (au./)^' are given 
simply in terms of those for £)^ M > j4 > ma P scaled by the 
appropriate factors of F^; A . 

The minimum-variance estimator &lm is then ob- 
tained by summing all of the individual estimators, for 
a given LM, with inverse-variance weighting (and tak- 
ing into account also the covariances). The minimum- 
variance estimator is thus 



LM, A, map 



V>1 AA' 
r 2 



and it has variance err, given by 

^ = E G «'E^r' he"')- 1 



V>1 AA' 



AA 1 



AA' 

(35) 
(36) 



In these expressions, the A A' sums are over 
{EB, BE, TB, BT, EE, TE, ET] for I ^ V and 
{EB, TB, EE, TE} for I = V, and the matrix inversion 
is in the AA' space. 



B. Variance for D 

We assume that the noise is isotropic and that it 
is not correlated with the data; that the temperature 
and polarization noises are not correlated; and that the 
temperature/polarization noises between different pixels 
are not correlated. If so, then the power spectra for 



the map are C\ 



A, map 



cfwf + C{ 



where here 



A= {TT, EE, BB}, as described in Ref. The TT, 

EE, and BB noise power spectra are 



a 



TT.noiso 



pix 



i or , 



s~iBB, noise s-iEE, noise / A / Af \ 2 
C t =C { = (47r/iV pix Jcrp , 



(37) 



where a\ and a 2 P are the temperature and noise variances 
in each pixel, and N p i x is the number of pixels in the map. 
The instrumental-noise contributions to the cross-power 
spectra are zero: 



EB: 



TB, 



TE: 



0. 



(38) 



Now we can calculate the desired variances for D, using 
Eqs. (J2J) and the power spectra Cf. For EB, for I ^ V, 
we have 



■,1V 



s-i B B ,ma,p ^1 E E ,m&p 



U BE,BE 

nil' B B ,map £jE E ,ma,p 



EB.EB 

II' 



= cp 



e — Cbe.eb — 0' 



and for I = V, 



pll ^BB.map^EE.map 

U BE,BE — °Z °Z 



(39) 



(40) 



The covariances for BT and TB are the same, with the 
replacements E— >T. There are also covariances between 
the TB and EB estimators. For I = V , they are 



For I ^ V, 



nil 

C BE,BT ■ 
nil' 

U BE,BT 
nil' 

U EB,TB 
ill' 



c, 



SB,map^fTE,map 



/-yBBjnap^TE^ap 



BB, map ^jTE, map 



nil nil 

^BE.TB — ^EB.BT 



I 

0. 



For EE, 



C 



w 

EE, EE 



(1 + <5h')Ci 



E E .map ^E E ,map 
°Z' 



For the TE case, for I ^ V , we have 



•IV 

TE,TE 
IV 

ET,ET 
II' 



^TT,map piE E , map 
°Z °Z' ) 

^TT^map^-f-Ei^map 
°Z' °Z i 



nil nil 

U TE,ET ~ C ET,TE 

and for 1 = 1', 



p-iT E , map piT E , map 
°Z °Z' ' 



C 



ET.ET 



£,TT, map £, EE, map ^(J TE 

,map 



(41) 



(42) 



(43) 



(44) 



(45) 



There are also TE-EE covariances. However, since TE 
and EE are almost always weaker probes of the rotation, 
we do not include these additional expressions here. 
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C. Quadratic Estimators for BE (or TE) only 

As an example (and for clarity), we can write down 
the expressions for the estimator and noise in the case 
where we use only information from the BE correlator to 
determine C(lm- In this case, the estimator is 



V>1 



jpL.BE ^.LM. BE, map 
IV U W 



c, 



BB, map ^jEE .map 



(B <-> E) 



(46) 



and the noise is given by 



V>1 



B B. map *~iEE. map 



(B <h> E) 



(47) 



The IV sums here are over values that satisfy I + I' + 
L =cven. The estimator and variance for TB are the 
same after the replacement E^T. 



We now derive the limit of high multipoles. 9 We start 
by using the approximation (2l| . 



i/jf, « cos 2ip u , 



I I' L 




(51) 



for the L + 1 + V =cven case. From Eqs. (J32J) and J5TJ, 
we have 



L,LV->oo IV f I V L 



-( ) cos 2 W . (52) 



From the relation of the spherical harmonics and Wigner 
3j symbols (see Appendix A), this gives, for large L, 



L,Z,Z'— >-oc 



IV 



2ttL 



(53) 



x / dhY m Y V0 Y L0 cos 2 2ip h 



Given that 



2 71 



n e — u m,Cb 

Z7T 



(54) 



IV. FLAT-SKY LIMIT AND COMPARISON 
WITH PREVIOUS WORK 

In this Section, we derive the flat-sky limit of the vari- 
ances for the rotation on the full sky and compare our 
results to those of Ref. [l5[ . We work out the EB case, 
where I + I' + L = even. The other three cases follow 
analogously. 

From Eqs. (|4"T|) and (f3~4")l . the variance is 



,r 2 = 4^(4)^(WiWr) 5 



V>1 



(CD 



{C EEy 



CB B ,map ^fE E ,map 
I U V 



Cr 
L i' 



BB, map ^jEE, map 



(48) 



The two terms in Eq. (|48[) are the same under the ex- 
change of I and V . Thus, after renaming the indices on 
one of the two terms, we get the sums over I < V and 
I > V, which covers the whole range of Vs a . We are left 
with 



K)- 2 = 4^4gJ,(4) 

w 

where we have defined 

Xt, = (W t Wv) 2 



s~tEE,niaps~<BB.map ' 



(49) 



(50) 



Eq. (|53p can be re-written using 



mm' M IV 



E E 

mm' M IV 



2nL 



dhY* m Y Vml Y LM 



(27T) 3 

We use relations from Ref. I21J| 



E* m ^ 



lm& 



6(L-(l-l'))= I * e'tf-**)* 



dft 



E YLYvm'Y LM < 



,'Af 

i(JW^Jx,-|-mV('— "ivO 



and replace the sum with the integral, 

E / <fy>j dtpv IV ^ ( { d 2 ld 2 V. 

j,f J J J 



(55) 



(56) 



W (57) 
ll'L 



(58) 



Note that, when we switch to integration, as shown further on in 
the text, the 1 = 1' term is of measure zero, and can be ignored. 



9 If L is large, then the triangle inequalities and the requirement 
for non-flat triangles ensures that / and V are also large. 
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From Eqs. (|55|) ([58)) and Eq. (f53]> . and after integrating 
over dtfL , we obtain the flat-sky limit for variance in the 
EB case, 



WMAP 



4 ' 



L.l.V 



x{wiWi-y 



(2tt)' 
{Cjf E f 



cos 2<pm 



(59) 



£iBB , map £iEE , map ' 



This can be shown to agree with the results of Ref. [15[ 
after combining their Eqs. (7) and (8) and results from 
their Table 1. 




V. NUMERICAL RESULTS 

We now present numerical results for the variances of 
the estimators for a position-dependent rotation, for dif- 
ferent instruments. The primordial power spectra, at the 
surface of last scatter, are obtained using WMAP-5 cos- 
mological parameters: Vl^h 2 = 0.02267, f2 c /i 2 = 0.1131, 
A = 0.726, n s = 0.960, r = 0.084, and a power spectrum 
normalized to WMAP5 [HI ]. 

We analyze three different experiments: (i) CMBPoPs 
(EPIC-2m) 150 GHz channel with resolution #f w h m = 5', 
taking the relevant parameters as given in Ref. [23[, 
a noise-equivalent temperature NET= 2.8/iK-yiec and 
the observation time t b s = 4 yr; (ii) Planck 143 GHz 
channel, with 0f w hm = 7.1', NET = 31/jK v / scc and 
t ohs = 1.2 yr; (iii) WMARwith 6> fwhm = 21', cr T = 
30/iK and o P = 42.6/^K 0. The NET parameters 
specified for Planck and CMBPol are related to the 
temperature/polarization pixel-noise variances through 
a^/N pix = / sky (NET) 2 /t obs , where / sky is the fractional 
sky coverage (here assumed to be unity). 

We evaluate the expressions for variances derived in 
the previous Section, for each of the three experiments, 
and show in Fig.[T]the numerical results for EB, EE, TB, 
and TE correlations. 

The first thing to notice about the variance levels in 
these three experiments is that both the lower noise and 
the higher resolution contribute to about two orders of 
magnitude improvement in the sensitivity to a rotation 
from WMAP to Planck and will lead to yet more than an- 
other order-of-magnitudc improvement in CMBPol. This 
is illustrated in more detail in Fig. [31 where the variances 
in the EB and TB estimators are compared for all three 
experiments. 

From Fig. [U we see that, at low multipoles (below L 
of about 200), the most sensitive estimators for all three 
experiments will be those derived from the EB and TB 
correlations. This comes about because of the absence of 
any TB or EB correlation under the assumption (justified 
largely by upper limits to the B mode amplitude) of no B 
modes at the surface of last scatter. We separately look 
at the TB-estimator variance for WMAP and the EB- 
cstimator variances for Planck and CMBPol, in Fig. [5J 



1 10 100 1000 

L 

Planck 

1000 r , -. 1 , , , , i i 




1 10 100 1000 

L 

CMBPol 




0.001 1 1 — 1 — 1 — 1 — 1 — 1 — 

1 10 100 1000 

L 

FIG. 1: Variances for the four rotation estimators for WMAP, 
Planck, and CMBPol are shown. We see that, at low mul- 
tipoles, the lowest-error estimator for Planck and CMBPol 
comes from the EB correlation, and for WMAP from the TB 
correlation. 
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FIG. 2: The variances of the estimators for rotation are shown 
in linear scale, for WMAP (variance from the TB correlator), 
Planck (from EB) and CMBPol (from EB). The range of the 
plots is chosen so that it covers the resolution domain of the 
instrument. The exponential rise at high L, due to the win- 
dow function, is visible. Note the difference in scale on the 
y-axis in all three panels. 



For WMAP, we find a TB-estimator variance of 8.3° 
at L = 0, which is consistent with the current constraints 
on a uniform rotation 10 [l(| EH (see Appendix C). The 
dipole and quadrupole components of the rotation have 
the same variances as the monopole [L = 0), since the 
variance remains fairly fiat out to L of about 100, in all 
three instruments. Above L ~ 400, the variance increases 
rapidly. This happens when the exponential part of the 
window functions dominates (due to the finite angular 
resolution of the instrument). Also, since the correlation 
angle for polarization is about 10 times smaller than that 
for temperature, the exponential tail in the EE-estimator 
case becomes prominent at higher multipoles than in the 
TE case. The variance from the EB-estimator at L = 
is 46.2° for WMAP, and thus not constraining. The TE- 
cstimator variance quickly drops below the EB-estimator 
variance (at L ~ 20) and below the TB-estimator vari- 
ance (at L ~ 150). However, the variance <jl at these L 
is so large (> 100) that the measurements are not at all 
constraining. Similar features are apparent in plots for 
the other two experiments. 

For Planck, the variances of the EB and TB estimators 
are more comparable, and the constraints to all rotation 
multipoles in the range from L of to about 300 come 
from the EB variance. At L = the variance is 14' 
and 24', for the EB- and TB-estimator variances, respec- 
tively. Planck can thus provide an order-of-magnitudc 
better sensitivity to the uniform rotation than the current 
WMAP sensitivity. For high multipoles, above L ~ 400 
or so, the TE-estimator variance becomes the smallest 
one. At L ~ 800, a rapid rise in all four variances is 
visible, due to the limitations in angular resolution and 
the correlation angle of the polarization. 

For CMBPol, the EB-estimator variance is the smallest 
in the whole range of multipoles from to 1000. At 
L = 0, the EB and TB values are, respectively, 2.2" and 
6.3" , which is better than Planck by more than an order 
of magnitude. Similarly to WMAP and Planck variances, 
we observe a rapid rise in the variance above L ~ 1000, 
corresponding to the resolution limitations and/or the 
polarization correlation length. 

In Fig. 2] we show the variance for the combined 
minimum-variance estimator, obtained from all four es- 
timators, with inverse-variance weighting. We have not 
included the covariance between the four estimators in 
this numerical calculation. However, this omission should 
make negligible difference for the WMAP and CMBPol 
curves, where the variance is determined primarily by TB 
and EB, respectively; the Planck curve may be increased, 
but only slightly. The run of this combined variance with 
L differs very little from the smallest /constraining vari- 
ance (that is TB in the case of WMAP and EB in the 



10 Note that we calculate the multipole moments of rotation, so a 
uniform rotation of angle a has aoo = v4?ra. Thus, ooo = 8.3° 
is equivalent to a = 2.3°. 
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FIG. 3: The best constraints to rotation, i.e. the variances for 
estimators from the EB and TB correlations, are compared for 
three instruments. Planck gives about 2 orders of magnitude 
improvement in comparison to WMAP, and CMBPol is one 
order of magnitude better than Planck. 
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FIG. 4: The combined variance, for all four estimators for 
rotation, is shown for WMAP, Planck and CMBPol. Note 
that, due to the inverse-weight summing, the smallest of the 
four variances dominates the shape of the curves. 



case of Planck and CMBPol), because that term domi- 
nates the sum. 

As a check, we repeated our calculations for the EB- 
corrclator using the flat-sky formulas from §IV, and com- 
pared those results to the full-sky. The two sets of vari- 
ances are in good agreement at high multipoles, where 
the flat sky makes a valid approximation (better than 
a fraction of a percent for L > 50). The discrepancy 
increases up to ~ 4% at lower L. Fig. \5\ compares the 
variances obtained from the full- and flat-sky treatments 



FIG. 5: Comparison between the full- and flat-sky expres- 
sions for the variance in the low-L regime. The discrepancy 
between the two is larger for even L because of the additional 
contributions to the estimators from 1 = 1'. 



at low L. Finally, we re-did the full-sky calculations using 
the same instrumental parameters as in Ref. (TsT ]; these 
are also in good agreement with our results at high L's 
(better than 7% for L > 50, and about 11% at L = 0). 



VI. CONCLUSIONS AND SUMMARY 

Exotic mechanisms, such as quintessence fields that 
couple to the pseudoscalar of electromagnetism, could 
rotate the linear polarization of the CMB. In this pa- 
per, we derive the complete set of minimum- variance es- 
timators for a position-dependent rotation of the CMB 
polarization, and thus provide a recipe, given a full-sky 
map, for measuring the rotation angle as a function of 
position on the sky. We also evaluate the variances with 
which the rotation-angle spherical-harmonic coefficients 
can be measured for WMAP, Planck, and CMBPol. Our 
results indicate that EB and TB correlations will provide 
more sensitive probes of the rotation angle than TE and 
EE correlations, and that EB becomes better, relative 
to TB, as the instrumental noise is reduced. We have 
checked that our results for the variances recover prior 
results, both analytically and numerically, in the flat-sky 
limit. As an additional check, Appendix C shows that 
our expressions for the rotation-angle variance recovers 
that expected for a uniform rotation. Appendix B shows 
that parity considerations can be used to distinguish the 
effects of rotation and of weak lensing on the CMB tem- 
perature/polarization map. 

Now that we have elucidated the all-sky formalism, the 
next step will be to work out algorithms for recovery of 
the rotation angle for a map with partial sky coverage. 
We anticipate that analogous techniques for determining 
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the cosmic-shear deflection angle may be adapted for this 
purpose. 
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The relation between spherical harmonics and Wigner 3j 
symbols is 

J dnY hmi (h)Yi 2 „ l2 (n)Yj 3m3 (n) 

/ (2Z 1 + 1)(2Z 2 + 1)(2Z 3 + 1) (70) 

x ( h h h \(h h h\ 

\ mi ?7i2 7713 Mo J 



APPENDIX A: Some Useful Formulas 

Some formulas used in §11 and §111 are included here. 
The metric tensor and its inverse on a unit 2-sphere (i.e. 
on the sky) is 
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1 
sin 2 



-1/10 

9 = ' -A- 



The Levi-Civita tensor on a unit 2-sphere is 



sin( 
-sinff 



An orthonormal basis on a unit 2-sphere is 









0- 




C ) 



(60) 



(61) 



(62) 



(63) 



Some useful properties of the Wigner 3j symbols and re- 
lated quantities include 



^i h h 
mi m>2 7773 



= (-1) 



h h h 
-mi — 7772 — 777-3 



h h h 
mi 777,2 77i3 



h h h 
7772 "73 mi 



h h h 
mi 77i2 1113 



= (-1) 



h+h+h 



I I 
777 — m 



mi + m,2 + 777.3 7^ 



(-1) 



h h h 
777,2 nil 7773 

l—m 
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h h h 
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(64) 
(65) 

(66) 
(67) 

(68) 
(69) 



APPENDIX B: Rotation vs. Weak Lensing 

The effects of rotation and weak lensing [2l], [24| on 
polarization are orthogonal and can thus be distinguished 
geometrically with a full-sky map. For example, if we 
start off with a pure E mode at the surface of last scatter, 
rotation induces a B mode, given by Eq. ([6]), where the 
only non-zero terms are those that satisfy L+l+l 2 =even. 
However, if we analyse the effect of weak lensing (see 
Ref. [l^l), a pure E mode polarization tensor changes by 



SP ab = (V c ^)(V c P ab ), 



(71) 



where (p is the projection of the gravitational potential 
along the line of sight. Thus, the B mode induced by 
weak lensing is 



SB hn = J dhSP ab (h)Y^ ah '(77) 

LM hrrt2 



dnY^\^ c Y (LM) ){^Y { l m2)ab ). 



(72) 



The parity of the spherical harmonic Ylm is (— 1) L ■ The 
parity of the E-mode term in the integral is (— l)' 2 , and 
the parity of the B-mode term is (— l) i+1 . The parity of 
the integrand in Eq. |jJ2J is then (— 1) z ' +f +' 2+1 . There- 
fore, the integral is non-vanishing only for terms that 
satisfy L + 1 + I2 =odd. 

We conclude that the rotation induces B modes that 
satisfy L + I + I2 =even, while for weak lensing we have 
L + l + h =odd [see Eqs. © and IpOjl]. Thus, the two 
effects can be entirely separated. 

Strictly speaking, this orthogonality between lensing 
and rotation occurs only at linear order in the rotation 
and lensing amplitudes in the limit that lensing and ro- 
tation are both small. If a lensed field is then rotated, 
and/or if a rotated field is then lensed, then the orthogo- 
nality will break down. However, this will occur only with 
an amplitude that is proportional to the product of the 
lensing and rotation amplitudes. We have here implicitly 
assumed this to be small and leave the full treatment of 
this higher-order effect for future work. 
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APPENDIX C: Uniform Rotation Cross-Check 

We can perform a cross-check of our formulas for vari- 
ances of the rotation, by analyzing only the L = term, 
where all the representation-theory coefficients can be 
readily evaluated. We do so for the EB case. From 
Eqs. (HHJ) and ([50)) (for I = I'), after evaluating the co- 
efficients for L = 0, we get the variance of the uniform- 
rotation estimator to be 



(°oo)~ 



J (2/ + l) 



^<££.map^(6B,map 



(73) 



We can see that this is indeed the right expression for 
uniform rotation, by noting that the B mode induced by 
small rotation, from a pure E mode, is given by Eq. (|3j). 
The induced EB power spectrum in that case is 



EB 



2aC, 



EE 
I ! 



(74) 



and estimators for the rotation can be expressed for each 
Ira pair as 



a 



pinap nmap 

2Cf E Wf ' 



(75) 



The variance of rotation is then calculated from all Im 
pairs as 

,2 



oc I 

-- 2 = EE 



1=0 ra=-l ( («) 



1 ,£ (21 + i) 



E E ,map ^jB B ,map 5 



(76) 

where the factor of (21 + 1) comes from the sum over m, 
because the terms in the sum have effectively only index 
I. For L = 0, Eq. (J7BJ) reduces to Eq. (73]), once the factor 
of s/4k by which a and aoo differ is taken into account. 



APPENDIX D: Contributions of Multipole pairs to 
Rotation Estimator 



In Fig. [SJ we show the summand in Eq. (|47p . versus I 
and I', for L = 20, 200, and 500. These terms are the 
weights with which each IV pair contributes to the EB 
estimator &lm for Planck (see Eq. (US])). Fig. [SJ there- 
fore illustrates which multipoles in the EE power spec- 
trum are expected to contribute most to the three chosen 
rotation multipoles. From the Figure, we see that the re- 
gion allowed by triangle inequalities grows with L, and 
that the peaks along I' correspond to the peaks in the 
EE power spectrum. There is also an overall exponential 
decay from the window functions, with a characteristic 
scale of I' ~ 1000, beyond which the distributions fall to 
zero. 

Overall, the EE multipoles that affect the estimate of 
the rotation at L = 20 and 200 seem to be predominantly 
those that correspond to the strongest peaks in the EE 
power spectrum, below the resolution limit of I' ~ 1000 




800 



FIG. 6: The summands in Eq. (|47[) are plotted versus I and 
I' . Plots for three different rotation multipole coefficients are 
shown: L = 20 (top panel), L — 200 (middle panel), and 
L = 500 (bottom panel). The plots suggest that the major 
contribution to any given rotation multipole coefficient comes 
from the 11' estimators that correspond to the strongest peaks 
in the EE power spectrum above the resolution limit of the 
instrument. 



(i.e. V of about 400, 700, and 1000; see Fig. [7J. The 
relative contributions are shown in the top two panels of 
Fig. [U The bottom panel of the same Figure shows the 
L = 500 case. There, we see that the largest contribution 
is local, i.e. it comes from non-flat triangles at multipoles 
below V ~ 500, and the secondary contribution comes 
from the EE power-spectrum peaks at V of about 700 
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FIG. 7: The EE power spectrum used for numerical estimates 
in §V is shown. The main peaks of the power spectrum are 
visible at L ~ 150, 400, 700, and 700. 
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